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ABSTRACT 



Context. The so-called shift parameter is related to the position of the first acoustic peak in the power spectrum of the temperature 

anisotropies of the cosmic microwave background (CMB) anisotropies. It is an often used quantity in simple tests of dark energy 

models. However, the shift parameter is not directly measurable from the cosmic microwave background, and its value is usually 

derived from the data assuming a spatially flat cosmology with dark matter and a cosmological constant. 

Aims. To evaluate the effectiveness of the shift parameter as a constraint on dark energy models. 

Metiiods. We discuss the potential pitfalls in using the shift parameter as a test of non-standard dark energy models. 

Results. By comparing to full CMB fits, we show that combining the shift parameter with the position of the first acoustic peak in the 

CMB power spectrum improves the accuracy of the test considerably. 
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1. Introduction 

Comparing cosmological models to current observational data 
can be cumbersome and computationally intensive. Large 
multi-dimensional parameter spaces cannot be probed by 
grid-based metliods but more sopliisticated approaches are 
required, for example Mon te Carlo Markov Chain s (MCMC) 
jGamerman & Lopes 20061 iLewis & Bridle 20021) . Of the 
most commonly used cosmological data sets available today, 
cosmic microwave background (CMB) and large-scale structure 
(LSS) observations in particular require computational effort 
in parameter estimation. Given computing time and patience, 
this is not a problem, at least in principle, when testing models 
in which the evolution of linear density perturbations is well 
understood and hence calculable. However, for several of the 
more imaginative dark energy models the situation is more 
complicated. In cases where the mod el is specified by an action, 
Uke the DGP model (iDvah et al. 2000 ), one should in principle 
be able to set up the equations for linear perturbations, but 
in practice this has turned out to be difficult and it is only 
recently that prog ress in this direction has b ee n made in this 
particular case (iKovama & Maartens 20061 iKoyama 2006t 



ISawicki et al. 20061: ISong et al. 20061) . And even though the 



equations are known, they may be so complicated to treat 
numerically as to make it practically impossible to explore 
the parameter space of the model in MCMC. In addition, 
there are a large number of dark energy models, based on 
phenomenological considerations, that lack the detail to allow 
one to proceed with well-defined calculations. Examples of 
such models include the various proposed modifications of the 
Friedmann equation, where the model is simply not specified 
well enou gh to allow the calcula t ion of the density perturb ations 
e.g. (preese & Lewis 20021: iGondolo & Freese 2003h . The 



(see 



justification of such models may be questioned, but the state 
of our understanding of dark energy argues for keeping an 
open mind. One would like to have some means of testing 
both groups of models, incorporating as much empirical in- 
formation as possible, but avoiding the need to calculate the 
behaviour of density perturbations. In practice, this means 
restricting the observational tests to those involving the age 
and distance scale, in particular the luminosity distance- 
redshift relationship as probed by supernovae of type la 
(SNIa). Important as the supernova data are, they are still 
not very restrictive if one allows for e.g. non-zero spatial 
curvature or a tim e -varying equation of s tat e for dark energy 



dRiess etal. 19981: iPerlmutter et al. 19991: iTonrv et al. 20oi. 



Barris et al. 20041; iRiess et al. 2004 lA stier et al. 200^ 



Wood-V asev et al. 2007 



Send offprint requests to: 0. Elgar0y 



Clocchiatti et al. 2()06t 

iMiknaitis et al. 20071: iDavis et al. 20071) ? 

To tighten up constraints on dark energy models, a common 
approach is therefore to include additional information about 
the distance scale from the CMB in the form of the so-called 
shift parameter that is related to the position of the first acous- 
tic peak in the power sp ectrum of the temperature anisotropies 
dEfstathiou & Bond 1999i) . Recently, after the baryon acoustic 
oscillations (B AO) where observed in the SDSS Luminous Red 
Galaxy sample (Eisenst ein et al. 2 005). it has also become com- 
mon to incl ude the information about the a ngular scale of the 
oscillations dPavis et al. 20071: IWright 200"7h . What one should 
bear in mind, however, is that these distance scales are not di- 
rectly measured quantities, but are derived from the observations 
by assuming a specific model, usually the flat ACDM model or 
a slight variation thereof. Care needs to be exercised when using 
these derived quantities to test more exotic dark energy models. 

Here we consider the shift parameter in more detail by com- 
paring its predictions to those obtained from full CMB fits for 
different types of cosmological models. We identify the limita- 
tions of using such a measure and advocate using a combination 
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of the acoustic peak scale along with the shift parameter as a 
more accurate probe of the CMB power spectrum. Such a com- 
bination is quick and easy to implement and should be included 
in tests of dark energy models where it is either cumbersome or 
unknown how to calculate the full CMB and matter power spec- 
trum. 

2. Theory 

The use of the shift parameter as a probe of dark energy is based 
on the observation that different models will have an almost 
identical CMB power spectra if all of the following criteria are 
satisfied (Efstathiou & Bond 1999 ): oj^ - D.ch^ and a>^ - Qb/i^ 
are equal, primordial fluctuation spectrum is unchanged, and the 
shift parameter, 

1/2 

K = ^sinnt(wyV), (1) 



1/2 



where sinn(x) 
tively, with 



{sin(x), X, sinh(ji;)) for k - 4-1,0,-1 respec- 
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is constant. In this original definition of the shift parameter, the 
universe is considered to be filled with matter (dark and bary- 
onic), w„ = ^tJt^i (J^hi curvature, u)^ - H^h^, cosmological 
constant wa - Oa/i^ and a dark energy component loq - Q.Qh^ 
with a constant equation of state w. The density parameter Q; is 
the ratio of the present-day density of component i to the density 
of a spatially flat universe, pc - 3//y/8;rG, and h is the dimen- 
sionless Hubble constant defined by Hq - lOO/z kms ' Mpc 
Integration is carried out from recombination, a^, until today 
a = 1. 

In a spatially flat universe (k = 0), the shift parameter re- 
duces to 

9? = Va: r ^ (3) 

Jo Eiz) 

where E{z) = H{z)/Ho and H{z) is the Hubble parame- 
ter. In this form the shift parameter has been used in a 
number of works ([L azkoz et al. 2006; Wang & Muk heriee 2006t 
lAmarzguioui et al. 2OO6:,Elgar0v & Multamaki 2005 1). 

The sound hor izon at recombination for three massless neu- 
trinos is given by ( lEfstathiou & Bond 1999h 

da 
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yj(a + a^^))(l + R(a)) 
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V^(«eq) 

where /?(a) = 30496wba andcieq = l/(24185wni) and the recom- 
bination^ redshift can be calculated by using the fitting formulae 
jHu & Sugivama 199@ : 

Zr = 1048(1 + 0.00124wb°'^^^)(l + giwg) (5) 
gi = 0.0783W|;"-2^V(1 +39.5w^™) 
g2 = 0.560/(1 + 21. It^i*^'). 

The OTth Doppler pe ak has the comoving wave number 
jHu & Sugivama 19961) mn - k,„rs{ai), and hence the location 
of the first peak in multipole space is approximately given by 



(6) 



where cIa - 'RluH^ is the angular diameter distance. Rewriting 
Eq. (|6]l, we have 



In 



1 + V/?(fleq) 



(7) 



The location of the first peak is hence a combination of the shift 
parameter and the size of the sound horizon at recombination, 
as is expected on physical grounds. Even though the relation be- 
tween 4 and is linear, the two parameters are not degenerate 
and in fact complement each other well in constraining models, 
as is shown later. 

The best fit value calculated from the WMAP-team provided 
MCMC chains for the shift parameter in the standard flat ACDM 
model is 



!R= 1.71 



0.03' 



(8) 



which is in good agreement with Wang & Mukherjee (2006). 
This result is practically equal for the ACDM and wCDM mod- 
els with or without dark energy perturbations. The acoustic peak 
position as measured by Eq. (|7]) calculated from the same data is 



4 = 303.6!|:^. 



(9) 



One should always bear in mind the conditions for the shift 
parameter to be applicable. If one wants to use the shift param- 
eter as a constraint on a dark energy model, then first of all 
the distribution of the shift parameter has to be derived from 
the CMB data. This cannot be done without assuming a model. 
Since one is only looking for a constraint on the expansion his- 
tory of the universe, it is easy to forget that one is also mak- 
ing assumptions about the primordial power spectrum of density 
fluctuations, since these form the basis for calculating the CMB 
anisotropics. In effect, one is therefore always making implicit 
assumptions about inflation, even though what one wants to test 
is the kinematics of the dark energy model under scrutiny. 

In order to demonstrate the significance of these underlying 
assumptions, we show the distributions of the shift parameter 
in figure [T] derived from MCMC chains for the ACDM model 
with four different primordial power spectra: the standard power- 
law version, power-law with running scalar spectral index, with 
tensor modes, and with both tensor modes and running scalar 
spectral index. The distributions for K are visibly different in the 
four cases. Therefore, whenever one uses the shift parameter one 
should be clear about the assumptions made in deriving its distri- 
bution from the CMB data. The acoustic scale 4, also shown in 
figure[Tl varies significantly less when changing the assumptions 
about the primordial power spectrum. 

In figure |2] we show the distributions for K and 4 when we 
assume a power-law primordial power spectrum, but make dif- 
ferent assumptions about the matter and energy content of the 
Universe. Here the distributions are less scattered, and again 
we note that the acoustic scale exhibits less variation than the 
shift parameter. As a demonstration of using the shift parameter 
blindly, we also plot the distribution of K and 4 when we allow 
for massive neutrinos. This is clearly false since the size of the 
sound horizon is now changed and hence the basis of using the 
shift parameter is no longer valid. This is important to take note 
of, because we know that neutrinos do have a mass that should 
always be included as a free parameter in cosmological parame- 
ter estimation, even though it is commonly neglected. 
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Fig. 1. The distribution of the shift parameter H and the acoustic 
scale 4 derived from MCMC chains with the WMAP data for 
four different types of primordial perturbations: power-law with 
no running scalar spectral index and no tensor modes (full line), 
power-law with tensor modes (dotted line), running spectral in- 
dex and no tensor modes (dashed line), and both tensor modes 
and running spectral index (dot-dashed line). 




Fig. 2. The distribution of the shift parameter H and the acous- 
tic scale 4 derived from MCMC chains with the WMAP data 
for four different models: ACDM (full line), dark energy with 
constant equation of state and dark energy perturbations (dotted 
line), open CDM (dashed line), and ACDM with massive neutri- 
nos (dot-dashed Une). 



3. Results 

Using the values of H and ^ calculated from the MCMC chains 
for the power-law wCDM model, we can compare the resulting 
confidence contours with those arising from doing the full CMB 
fit. 



3.1. w CDM model 

In figure[3]we show the 68%, 95% and 99% confidence levels 
arising from using the shift parameter and the acoustic peak posi- 
tion for the flat wCDM model. In calculating the confidence lim- 
its using the acoustic peak position, we have chosen a flat prior 
h - 0.73±0.03 and marginalized over h. We have kept the baryon 
density fixed at the WMAP value (jl>\, - 0.0223 which we use 
throughout the paper unless otherwise stated. Comparing this 




Fig. 3. The 68%, 95% and 99% confidence levels calculated 
from the shift parameter and acoustic peak position 



with the probability density plot from the full MCMC chains, 
e.g., for the wCDM model with no dark energy perturbation, fig- 
urelH we see that the shift parameter gives a good approximation 
to the CMB data while the acoustic peak position does not. Note 
that even though is fixed, H and /"a contours have fundamen- 
tally different shapes, demonstrating the importance of the sound 
horizon size in calculating l^. 




Fig. 4. Normalized probability density function from the WMAP 
3-year data 



The fact that the shift parameter approximates the full CMB 
contours so well for the wCDM model is not surprising since 
the value of the shift parameter has been calculated from a chain 
that assumes that the cosmology is of the wCDM type. In other 
words, we first assume a model and then calculate chains that 
best fit the data from which we derive a quantity. Following the 
same prescription one can in fact construct other quantities that 
also well approximate the full CMB contours, but are not physi- 
cally motivated. 

3.2. Role of the Hubble parameter 

In the WMAP chains, the Hubble parameter typically has a fairly 
large prior, 0.5 < h < 1 .0. If the value of h is constrained by other 
observations, e.g., the Hubble Key Project reports h - 0.72+0.08 
(Freedmanet al. 2001), the contours in the (Q^, w)-plane look 
quite different. In figure |5] we show the normalized probability 
density for the wCDM model from the WMAP provided chains 
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(no dark energy perturbations) with cuts on h along with the con- 
fidence contours calculated by using {.^ with and without the shift 
parameter. The shift parameter is independent of h and hence the 
confidence contours are unchanged and shown in figure[3] From 




Fig. 5. Comparison of different methods: probability density 
from the wCDM MCMC chains (top left), confidence contours 
calculated using {.^ (top right) and the combined 4 contours 
(bottom). 



the figure one can conclude that the combination of K and ap- 
pears to be a good approximation to the full CMB data when h 
is constrained by independent observations. This is further sup- 
ported when we consider the confidence contours arising from 
using the CR, 4) combination with h - 0.5 -1.0 (flat prior), 
shown in figure|6l Comparing to figure|4]we see that also in this 




effect on the results. In particular, when compared to the effect 
of changing the Hubble parameter, the significance of varying 
the baryon density within the BBN limits is negligible. 



4. Non-standard cosmologies 

The shift parameter is particularly useful as a quick measure of 
how a given cosmological model fits the CMB data. In order to 
assess the validity of this approach, we compare here the param- 
eter constraints arising from the shift parameter and from doing 
the full CMB fit on a non-standard model, namely on a general 
Friedmann equation. Such a model is generalization of the stan- 
dard Friedmann equation and as such serves as a useful generic 
no-standard model. We also consider the use of the acoustic peak 
position as an useful approximation to the full CMB data fit. 



4.1. Modified Friedmann equation 

Modified Friedmann equations arise, e.g., in alternative theo- 
ries of gravity. As an example, in the well known DGP model 
dPvali et al. 20001) . the Friedmann equation on the brane is of 
the form 



H 8kG 



-Pn 



(10) 



where pm is the matter density on the brane and rc is the cross- 
over scale at which gravity starts to feel the effects of the fifth 
dimension. 

Here we will consider modifications to the Friedmann 
equation with no spatial curvature in the spmt of 
jDvali & Turner 20031; lEigar0v & Multamaki 2005h . A general- 
ized Friedmann equation can be written as 



f{H,H,)^Hl^lM+z? 



(11) 



where instead of modifying the matter content we consider mod- 
ifications of gravity by having an arbitrary function /. The crit- 
ical scale, //c, is close to the present Hubble parameter, Hq, and 
determines when modifications from the standard Friedmann 
equation start to have an effect. At early times, when H » H^, 
we know from BBN constraints that /(//) ~ H^. Keeping this 
mind and expanding in terms of He I H gives 



Fig. 6. The 68%, 95% and 99% confidence levels calculated 
from the combination of the shift parameter and 4, A = 0.5 - 1 .0 



case, the combination of and 4 is a reasonable approximation 
to the results obtainable from doing the full CMB fit, figure|4] 



3.3. Role ofbaryons 

In the previous calculations, we have kept the baryon density 
fixed at the WMAP 3-year value, Wb = 0.0223. This value is de- 
rived assuming the ACDM model and hence when using the or 
4 to study other cosmologies, one should be somewhat cautious 
when using this value. A more robust, with respect to changing 
cosmology, measure of Wb comes from Big Bang nucleosynthe- 
sis (BBN), which gives 0.017 < Wb < 0.024 (95% confidence 
level) with three massless neutrinos. We find that changing the 
baryon density within the 95% limits from BBN has only a small 



1 + 



(12) 



from which it is clear that non-standard effects only start to have 
an effect at late times when H ~ H^. Expanding the sum gives 



1 +c 



H \H I 



^H^QM +zy. 



(13) 



In this form, one can interpret the cosmological constant as a 
second order correction to the Friedmann equation while the first 
order correction corresponds to the DGP model. Generally, the 
nth order correction for a flat universe is hence 



=Q^(l+z)^ + (l-Q™)(^)" 



(14) 



The leading correction to the Friedmann equation was previ- 
ously studied using current CMB, SNIa and LSS data by Elgar0y 
and Multamaki (2005). 
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4.2. Confidence contours 

The CMB spectrum arising within the context of the modi- 
fied Friedmann equation, Eq. (fT4l i. can be straightforwardly cal- 
culated by using, e.g., CMBFAST ( Seljak & Zaldarriaga 1996) 
(see Elgar0y & Multamaki (2005) for a detailed description) 
and fitted to the WMAP 3 year data. The resulting confidence 
contours are shown in figure Q along with contours from using 
the shift parameter and the acoustic peak position. In figure [T] 
we also show the contours obtained by combining the two ap- 
proximate methods. Again, we find that the combination of the 




Fig. 7. The 68%, 95% and 99% confidence levels calculated 
from the shift parameter (top left), acoustic peak position (top 
right), the combined K and (bottom left) contours and from 
the full WMAP 3yr TT data (bottom right) {h = 0.73 ± 0.03, flat 
prior) 



acoustic peak along with the shift parameter gives a good ap- 
proximation to the fuU CMB fit. 

5. Conclusions 

In this work we have reconsidered the use of the shift parame- 
ter as a quick and easily implementable probe of dark energy. 
We find that while it gives an excellent measure of the CMB 
spectrum for the ACDM model, with Wb and llIc fixed, caution 
should be exercised when using the shift parameter to compare 
and constrain non-standard cosmological models. 

Even when using to constrain more standard type models, 
such as the wCDM model, careful consideration should be given 
to the value of the Hubble parameter, h. Since the shift parameter 
is independent of h, but the full CMB spectrum fit is most def- 
initely not, the shift parameter can be misleading when applied 
blindly. 

The shift parameter is a geometrical measure as it measures 
the size of apparent sound horizon at recombination. Keeping 
the sound horizon size fixed, different cosmological models lead 
to different background expansion and hence the shift parameter 
can be used to compare and constrain different models. However, 
also the sound horizon size changes when varying cosmological 
parameters, most notably changing the matter density, Qm and 
the Hubble parameter, h. In addition, massive neutrinos will also 
have an effect. Hence, in general the shift parameter will not be 
an accurate substitute for the CMB data and may in principle 
give misleading results when used to constrain non-standard re- 
sults. 



In addition, the value of the shift parameter used to con- 
strain different dark energy models is derived by first assuming 
the ACDM model, fitting the model to the data and then cal- 
culating the value of K. Again, for a general model, using the 
value obtained in this manner is questionable since for a differ- 
ent model one may expect the shift parameter to be different, 
while the CMB spectrum can fit well with observations. 

In order to enhance the effectiveness of using the shift pa- 
rameter as a cosmological tool, we have considered adding infor- 
mation from the location of the first CMB peak, Combining 
these two easily calculable observables, allows one to encom- 
pass information from both the size of the sound horizon at re- 
combination and the angular diameter distance to it. As such, it 
more effectively constrains the allowed parameter space, includ- 
ing the Hubble parameter that is not fixed when using only the 
shift parameter. A possible caveat is, again, the fact that the nu- 
merical value of both of these parameters is calculated within 
the ACDM framework, but by comparing to different models, 
we see that the combination proves to be a good and efficient 
probe of non-standard cosmologies. 
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